Mon. Not. R. Astron. Soc. 000, [THU] (0000) Printed 17 December 2008 (MN I^TeX style file v2. 2) 



Light element abundances in carbon-enhanced metal-poor 
stars 



Richard J. Stancliffe* 

Centre for Stellar and Planetary Astrophysics, Monash University, PO Box 28M, Clayton VIC 3800, Australia 



Accepted 0000 December 00. Received 0000 December 00; in original form 0000 October 00 



ABSTRACT 

We model the evolution of the abundances of light elements in carbon-enhanced metal- 
poor (CEMP) stars, under the assumption that such stars are formed by mass transfer 
in a binary system. We have modelled the accretion of material ejected by an asymp- 
totic giant branch star on to the surface of a companion star. We then examine three 
different scenarios: one in which the material is mixed only by convective processes, 
one in which thermohaline mixing is present and a third in which both thermohaline 
mixing and gravitational settling are taken in to account. The results of these runs 
are compared to light element abundance measurements in CEMP stars (primarily 
CEMP-s stars, which are rich in s-processes elements and likely to have formed by 
mass transfer from an AGB star), focusing on the elements Li, F, Na and Mg. None of 
the elements is able to provide a conclusive picture of the extent of mixing of accreted 
material. We confirm that lithium can only be preserved if little mixing takes place. 
The bulk of the sodium observations suggest that accreted material is effectively mixed 
but there are also several highly Na and Mg-rich objects that can only be explained 
if the accreted material is unmixed. We suggest that the available sodium data may 
hint that extra mixing is taking place on the giant branch, though we caution that 
the data is sparse. 

Key words: stars: evolution, stars: AGB and post-AGB, stars: carbon, stars: Popu- 
lation II, binaries: general 



1 INTRODUCTION 

Carbon-enhanced metal-po or (CEMP) sta rs are defined as 
stars with [C/FeQ>-l-1.0 (|Beers fc Christl icb 2005), with 
[Fe/H] < —2 in most cases. These o bjects appear with in - 
creasing frequency at low metallicity l|Lucatello et al.l|2006l ) . 
The study of CEMP stars is being used to probe con- 
ditions in the early universe. For example, CEMP stars 
have been used to infer the initial ma ss function in the 
early Galaxy (e.g. iLucatello et al.ll2005l ). Chemical abun- 
dance studies have revealed that the majority of the CEMP 
stars are rich in s-process elements like barium l|Aoki et al.l 
I2OO3I ). forming the so-called CEMP-s group. Recent survey 
work has detected radial velocity variation in around 68% 
of these CEMP-s stars a nd this is consistent w ith them all 
being in binary systems (jLucatello et al.ll2005l ). 

Binary systems provide a natural explanation for these 
objects, which are of too low a luminosity to have been 



able to produce their own carbon. The primarjj^ of the 
system was an asymptotic giant branch (AGB) star which 
became carbon-rich through the action of third dredge-up 
(the deepening of the convective envelope into regions of 
the star where material h as undergone nuclear burning, see 
e.g. Iben & Rcnzini 1983) and transferred material on to the 
low-mass secondary (most likely via a stellar wind). The pri- 
mary became a white dwarf and has long since faded from 
view, with the carbon-rich secondary now being the only 
visible component of the system. 

It has commonly been assumed that accreted material 
remains unmixed on the surface of the recipient star during 
the main sequence and only becomes mixe d with the stel- 
lar interior during first dredge-up. However. IStancliffe et al] 
(2007) pointed out that this picture neglects thermoha- 
line mixing. Thermohaline mixing is the process that oc- 
curs when the mean molecular weight of the stellar gas in- 
creases toward the surface. A gas element displaced down- 
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^ This is the initially more massive star in the system; the sec- 
ondary is the initially less massive star. 
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wards and compressed will be hotter than its surround- 
ings. It will therefore lose heat, become denser and continue 
to sink. This leads to mixing on thermal ti mescales until 
the molecular weight difference is eliminated (jUlrichlllQTj : 
iKippenhahn et al]|l980l) . Accreted material has undergone 
nuclear processing in the interior of the companion AGB star 
and hence has a higher mean molecular weight than the pris- 
tine material of the companion. The accreted AGB material 
should therefore become mixed with the pristine material of 
the secondary b ecause of the action of thermohaline mixing. 
In their models, IStancliffe et al.l |20o3) showed that thermo- 
haline mixing could, under optimal circumstances, lead to 
the accreted material being mixed with nearly 90 per cent 
of the secondary. 

The extent to which accreted material is mixed into 
the accretor during the main sequence has been questioned. 
Using a sample of barium-rich CEMP stars (i.e. a set of 
stars where we expect t he AGB mass transfer scenario to 
hold), lAoki et al.l l|2008l ) showed that the distribution of 
[C/H] values in turn-off stars (i.e. those stars that have 
reached the end of their main-sequence lives, are still of 
low-luminosity and have yet to become giants) was differ- 
ent from that in giants suggesting that significant mixing 
onl y happened at first dredge-up . A sim ilar point was made 
by Denissenkov fc Pinsonneauhl l|2008bl ) using the data of 
iLucatello et al. ( 200&) . These authors showed that the abun- 
dance patterns of turn-off stars and giants were consistent 
with coming from the same distribution if first dredge-up 
resulted in the [C/H] value (and also the [N/H] value) being 
reduced by around 0.4 dex. This drop in both C and N is the 
result of the accreted material being mixed with pristine stel- 
lar material during the deepening of the stellar envelope at 
first dredge-up. They find that this result is consistent with 
having an accreted layer of material mixed to an average 
depth of about 0.2 M0 (or alternatively having an accreted 
layer of 0.2 M0 that remai ns unmixed until first dredge- 
up). [Thompson et al] (|2008l ) also questioned the efficiency 
of mixing in the binary system CS 22964-161. This system 
consists of two turn-off stars which are both carbon-rich. The 
system is also found to be lithium-rich, with e(Li) — +2.0^ 
As lithium is such a fragile element and is easily destroyed 
at temperatures above around 10® K, any extensive mixing 
of accreted mate r ial w ill result in lithium being depleted. 
[Thompson et al. I (|2008') suggest that gravitational settling 
could in principle inhibit the action of thermohaline mixing 
as helium diffusing away from the stellar surface would cre- 
ate a stabilising mean molecular weight gradient (a so-called 
'/x-barrier'). 

This idea was tested bv IStanchffe fc Glebbee^ (|2008l ). 
who included the physics of both thermohaline mixing and 
gravitational settling in their models of stars accreting ma- 
terial from a putative AGB donor. Their grid of mod- 
els over a range of accreted masses (0.001-0.1 Mq) from 
donors of diff'erent initial mass (1-3.5 M0) showed that grav- 
itational settling could only seriously inhibit thermohaline 
mixing if a small quantity o f material had been accreted. 
IStancliffe fc Glebbeekl (|2008l ) compared their predicted sur- 
face abundances of carbon and nitrogen to those of observed 



^ e(Li) = logj^Q ^.^1^^4-12 where TV is the number abundance of 
the element. 



CEMP stars, concluding that none of their model sets did a 
good job of reproducing the observed abundances patterns. 
The aim of this work is to extend the secondary models of 
SG08 to a greater range of isotopes (rather than just the G 
and N looked at in that paper) in the hope that these ad- 
ditional isotopes may help to illuminate the processes that 
happen in the AGB mass transfer formation scenario. 

Recently, there has been considerable progress in mea- 
suring the abund ances of light elem ents in carbon-enhanced 
metal-poor stars. lAoki et "all (|2008l ) have provided a sample 
of over 70 barium-rich CEMP stars that have other abun- 
dance determinations, including for sodium and magnesium. 
The element fluorine h as also been m easured for the first 
time in a CEMP star jSchuler et al.l [2007, 1. Each element 
that we have measurements for potentially gives us another 
way to explore both the nature of the donor AGB stars and 
what happens to the material that is accreted on to the sec- 
ondary star. This paper looks at what the light elements in 
CEMP stars can tell us about these topics. The AGB mass 
transfer scenario is likely to apply to the majority of CEMP 
stars (specifically those with s-process enrichments which 
are classified as CEMP-s stars), but not all CEMP stars 
form this way. In addition, the elements discussed in this 
paper do not necessarily come from AGB stars (for example, 
Na and Mg in the Galaxy mostly come from supernovae). 
The reader should bear these points in mind throughout. 



2 THE STELLAR EVOLUTION CODE 

Calculations in this work have been carried out using a mod- 
ified versio n of the STARS ste llar evolution code originally de- 
velo ped bvlEggletonl (|l97lf ) and updated by many authors 
fe.g. lPols et al. 19951 ). The code solves the equations of stel- 
lar structure and chemical evolution in a fully simultaneous 
manner, iterating on all variables at the same time in order 
to converge a model fsee IStancliffdlioO^ . for a detailed dis- 
cussion). The version used he r e incl udes the nucleo s ynthe - 
sis routines of IStancliffe et al.1 (120051 ) and IStanch^ l|2005h . 
which follow the nucleosynthesis of 40 isotopes from D to 
■^■^S and important iron group elements. Th e code uses the 
opacity routines of lEldridge fc ToutI (|2004|). which employ 
interpolation in the OPAL tables (jlglesias fc Rogerslll996h 
and which account for the variation in opacity as the C and 
O content of the material varies. 

To produce our model s, we follow the procedure of 
IStancliffe fc Glebbeekl (|2008l ). hereinafter SG08. We take the 
average composition of the ejecta (selected abundances are 
shown in Table [1} from the 1.5 Mq model of SG08 and ac- 
crete it on to our secondaries at a rate of 10~® Mq yr~^. The 



SG08 AGB mode ls have been evolved with the same code 
used here (see e.g. IStancliffe et al ] |2004 for details), so these 
models remain fully self-consistent. In each case, we accrete 
enough material to make the final mass of the star equal to 
0.8 M0, which is the appropriate turn-off mass for the halo. 
The accretion event takes place at a time appropriate to the 
age at which the primary enters its superwind phase. We 
use initial secondary masses of 0.7, 0.75, 0.78, 0.79, 0.795, 
0.798 and 0.799 M© and the metallicity is Z = 10"" (cor- 
responding to [Fe/H] = —2.3). Three model sets have been 
evolved: a 'standard' model which includes only canonical 
mixing processes (i.e. convection), the second includes ther- 
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mohaline mixing of the accreted material and the third se- 
quence includes both thermohaline mixing and gravitational 
settling. Mass loss has not been taken into account in the 
evolution of these models. 



3 RESULTS 

One caveat of the models must be borne in mind through- 
out. The AGB models from which the nucleosynthesis is 
taken and the detailed models of the secondaries presented 
here have been calculated using solar-scaled abundances and 
opacities. The former should only be important for sec- 
ondary species (i.e. those that derive from metals already 
present in the star at the time of formation). Most of the 
species that we are concerned with are primary and have 
been manufactured from hydrogen and helium within the 
star. The yields of these primary isotopes should be inde- 
pendent of the choice of initial abundance. As for the opac- 
ities, it would be more appropriate to use a set of tables 
that accounts for the enhancements of the Q-elements ex- 
pected in low-metallicity objects. However, an a-enhanced 
set of tables with variable C and O composition as used by 
the code is not currently available. We would expect that 
an a-enhanced mixture (with increased abundances of ^''O, 
^'^Ne and ^^Mg) would affect the yields of ^^Ne, ^^Na and 
25,26-|y|-g ^jjg production pathways for these elements are 
quite complex and interdependent, it is not immediately ob- 
vious how they would change. Computation of such models 
merits future study and is beyond the scope of this work. 



3.1 Lithium 

Figure [1] shows the evolution of the lithium abundance as 
a function of luminosity for a selection of the models. The 
standard models show the least depletion of lithium because 
the accreted material remains in layers that are too cool 
to destroy this fragile element. The only reduction comes 
at first dredge-up when the deepening convective envelope 
mixes the accreted material with the interior of the star. 
Models including only thermohaline mixing show the great- 
est degree of depletion, with the reduction being greatest 
when the most mass is accreted. This is because material is 
mixed deep into the star, dragging lithium that would nor- 
mally (in the absence of any mixing) remain near the stellar 
surface to depths where it can be destroyed. In fact, the pres- 
ence of Li in CEMP stars could be the best test to determine 
whether mixing has occurred after the accretion process. A 
high lithium abundance precludes the occurrence of efficient 
thermohaline mixing alone (however, see below). The addi- 
tion of gravitational settling alleviates the destruction of Li 
somewhat, with a drop in the degree of reduction. In the 
case that only a small quantity of material is accreted, the 
abundance can be preserved at close to the accreted value 
(though gravitational settling causes this value to drop over 
the main sequence). 

The initial lithium abundance used in both the AGB 
models and the secondaries are most likely too low. We 
should expect the i nitial value to reflec t the Spite plateau 
value of e(Li) « 2.3 (jSpite fc Spite|[l98^ ') . However, this may 
not matter for the AGB models, as lithium should easily be 
destroyed in the course of a star's evolution. In fact, in the 
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Figure 1. The evolution of the lithium abundances as a func- 
tion of luminosity for each of the model sequences. The standard 
models are displayed with solid lines and those with only thermo- 
haline mixing are displayed with dashed lines. The models with 
both thermohaline mixing and gravitational settling are shown by 
dotted lines. From top to bottom, the panels represent the case 
of accreting 0.002, 0.01 and 0.1 Mq of material respectively. In 
each case, material has been accreted from a 1.5 Mq companion 
and the final mass of the star is 0.8 Mq. 



1.5 Mq model, the lithium abundances drops by over two or- 
ders of magnitude on the first giant branch. For the 1.5 Mq 
accreted material, we would not expect the AGB star to have 
produced any significant amounts of Li during its TP phase 
(although see Abia &^Isern 1997 and Uttcnthalcr ct al. 20(53 
who report the detection of Li in low-mass Galactic AGB 
stars). Would an enhanced Li abundance in the secondary 
make a difference? If mixing is extensive, then Li will be de- 
pleted and an enhanced 'resevoir' in the surface layers will 
not matter. If the mixing of accreted material is shallow then 
an enhanced Li abundance in the secondary would lead to a 
higher abundance after mixing. 

Large lithium abundances have been measured in some 
CEMP stars. For example, CS 22964-161 is a binary CEMP 
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1.5 


0.690 


1.22(-3) 


2.78(-5) 


1.79(-11) 


8.83(-7) 


5.09(-6) 


3.81(-6) 


4.00(-6) 


2.23(-6) 


5.74(-6) 


2 


0.678 


1.62{-2) 


2.61{-5) 


1.18(-11) 


7.95(-7) 


1.66(-5) 


8.56(-6) 


2.45(-5) 


2.17{-5) 


5.61(-6) 


3.5 


0.679 


2.09{-3) 


6.57{-3) 


6.45(-13) 


2.14(-9) 


5.89(-5) 


5.25(-6) 


2.72{-5) 


7.55{-5) 


5.70(-6) 


Secondary 


0.749 


1.73{-5) 


5.30{-6) 


4.68(-ll) 


2.03(-9) 


1.67{-7) 


2.57(-6) 


3.83{-7) 


3.88{-7) 


5.85(-6) 



Table 1. The abundances of selected isotopes in the AGB ejecta used in these simulations along with the initial abundances used in the 
secondary models. The abundances are in the format n{m) = n X 10™. 



system which shows a lithium abundance of e(Li) = 2.09 
iThompson et al.llioOSl 'l. The AGB models do not have suf- 
ficiently iiigii lithium abundance to match this value, even 
before any mixing of the accreted material takes place 
(whether this be thermohaline mixing after accretion or con- 
vective mixing during first dredge- up). While it is possible 
for the more massive AGB stars to produce lithium via the 
Cameron-Fowler mechanism, such a process does not oper- 
ate in the low-mass AGB stars. We would have to suppose 
the existence of some additional mixing mechanism, active 
during the AGB, which would cause circulation of mate- 
rial below the base of the convective envelope down to tem- 
peratures where lithium production can take place. Such a 
supposition may be reasonable in light o f the existence of Li- 
rich low-mass AGB sta rs in the Galaxy (|Abia fc Isern|[l997l : 
lUttenthaler et al.ll2007l 'l. 

Attention should also be drawn to the Li-rich star 
H KII 17435-00532 . This star has e(Li) = 2.1 as reported 
bv lRoederer et al] l|2008h . The exact evolutionary status of 
this object is unclear, but it is clear that it has evolved at 
least as far as the red giant branch - yet it is still Li-rich! 
Roederer et al. suggest that it cannot have acquired its Li 
from a companion: other elements show evidence for a large 
degree of dilution. This would require the accreted mate- 
rial to have an unreasonably high Li abundance. They dis- 
cuss potential sources for Li production, but these all occur 
at late stages in the stars evolution. Unlike HKII 17435- 
00532, CS 22964-161 is still at a relatively early stage in 
its evolution and could not have undergone the proposed 
Li-enrichment events. We must therefore seek a different ex- 
planation for its Li-enrichment. 

Recent observations of stars in the metal-poor 
([Fe/H] = -2) gl obular cluster NGC 6397 (|Korn et al.ll2007l : 
iLind et al.ll2008l ) suggest that models that include gravita- 
tional settling are able to reproduce the abundance trends 
of certain element s if an ad ditional turbulent mixing process 
is included (Ric hard et al.|[2005 ). We have therefore run a 
model sequence ad opting the following prescription (from 
iRichard et aPbOOSl ) for such an extra mixing process: 



D = 400Dho{To) 



Pin) 



(1) 



where D is the turbulent diffusion coefficient, Z)He(Tb) is 
the atomic diffusion co efficient for hel i um at a temperature 
of To and p is density. I Richard et al.l (|2005 h conclude that 
a value of To = 10^ K can reproduce the Spite Plateau in 
Population II stars. We have therefore adopted this value. 
Note that there is (as yet) no physical cause for such mixing 
and the prescription remains an ad hoc one. 

We have recomputed the 0.1 M© case with this addi- 




Figure 2. The evolution of the lithium abundance as a function 
of luminosity when 0.1 Mq of material is accreted. The solid line 
shows the case where thermohaline mixing and gravitational set- 
tling are included. The dashed line has the additional turbulent 
mixing of lRichard et al" 1 l l2005D included. 



tional mixing mechanism included (in addition to thermo- 
haline mixing and gravitational settling). The results are 
shown in Figure [S] The addition of extra mixing reduces the 
height of the /i-barrier and hence allows somewhat more effi- 
cient thermohaline mixing. The degree of Li depletion across 
the main sequence is significantly reduced because the ad- 
ditional mixing inhibits gravitational settling by stirring up 
the surface regions. The extra mixing does not substantially 
change the Li abundance however, as extensive mixing can 
still take place and it is this that causes the destruction of 
Li. 



3.2 Fluorine 

The evolution of the fluorine abundance for a selection of 
the models is showing in Figure [3] Without any additional 
mixing mechanism, the secondary retains the surface abun- 
dance of the AGB ejecta which is [F/Fe] ~ 2.3. It then suf- 
fers a dilution at first dredge-up which, depending on the 
mass of material accreted, can cause [F/Fe] to drop by be- 
tween 0.5-2 dex. Models including thermohaline mixing have 
an immediate dilution of fluorine of around this level and 
may also suffer a further dilution at first dredge-up, depend- 
ing on whether the accreted material is mixed to a depth 
greater than the depth of the convective envelope during 
first dredge-up. In the case that both thermohaline mixing 
and gravitational settling are included, the fiuorine abun- 
dance suffers less of a reduction on the main sequence, but 
only if a small quantity of material is accreted. A surface 
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Figure 3. The evolution of [F/Fe] as a function of luminosity 
for each of the model sequences. The standard models are dis- 
played with solid lines and those with only thermohaline mixing 
are displayed with dashed lines. The models with both thermoha- 
line mixing and gravitational settling are shown by dotted lines. 
From top to bottom, the panels represent the case of accreting 
0.002, 0.01 and 0.1 Mq of material respectively. In each case, ma- 
terial has been accreted from a 1.5 Mq companion and the final 
mass of the star is 0.8 Mq. 

abundance of [F/Fe] of around 2 or less is predicted. This 
value declines slightly over the main sequence due to the 
effects of gravitational settling. 

Fluorine has only been measured in one CEMP star 
to date: HE 1305-1-0132. This object has an extremely su- 
perso lar F abundance of [F/Fe] = -|- 2.90 (jSchuler et al.l 
l2007l 'l. The metallicity is [Fe/H] = -2.5 and the s tar does 
show barium lines in its spectrum (|Goswamil |2005^ and so 
we should expect the model presented here to be appropri- 
ate. HE 1305+0132 is a red giant with a surface gravity of 
logio 9 ~ 0-8 S'lid hence we should compare this to the mod- 
els once they have passed through first dredge-up. Unfortu- 
nately, we are unable to distinguish between the three model 



sequences once the star has passed through first dredge-up 
as they all give similar [F/Fe] values. However, it is difficult 
to reconcile the models with the observations. There is insuf- 
ficient fluorine in the AGB ejecta to start with as the models 
do not reach high enough [F/Fe] on the main sequence, let 
alone after dilution of material has taken place either via 
thermohaline mixing or during first dredge-up. The case is 
worse if less material is accreted. If 0.1 Mq of material is ac- 
creted the models are about 1 dex below the observed [F/Fe] 
of HE 1305+0132, while if only 0.001 Mq are accreted the 
discrepancy is over 2 dex. However, the [C/Fe] and [N/Fe] 
values are quite well reproduced by a model that has ac- 
creted 0.1 Mq of material from a 1.5 Mq companion and 
mixed this material during the main sequence (this is neces- 
sary to elevate the nitroge n abundance, which h appens via 
CN-cycling as describe bv lStanchffe et al.1 l|2007l ).'). On the 
giant branch, this model has [C/Fe] = 2.0 and [N/Fe] = 1.5, 
which compare very favo urably with the mea sured values of 

2.2 and 1.6 respectively (jSchuler et aLlbOOTl ). 

The only way out of this dilemma is to suggest that the 
AGB models are underabundant in fiuorine. There could be 
two causes for this: the model has a companion of a too low 
a mass to produce enough fiuorine, or there is an intrinsic 
problem with the pr oduction of fiuorine i n the AGB models. 
On the former point. Iliugaro et al.l (|2008l ) have shown that a 
2.25 Mq star is likely to produce the maximum fluorine yield. 
On the latter point, the nucleosynthesis of fluorine takes 
place via an involved chain of reactions (see iLugaro et al.l 
|2004 for details). There are many uncertainties associated 
with the production and destr uction of '^^F (see the dis- 
cussion in 'Karaka s et al1l2008l ). On top of this, there are 
mod el uncertainties such as the efficiency of third dredge-up 
(e.g. iFrost fc Lattanzidll996l : iHerwid [2OO0I : IStanchffe et"al] 
20051. among m any others) and the mass loss rates used (e.g. 
Stancliffe fc Jeff crv 20070. As an indication of such eff ects, 
we note that the models of lKarakas fc Lattanziol l|2007l ) are 
more abundant in ^^F than the SG08 models used here, by 
roughly an order of magnitude at 2 Mq . This may allow the 
reconciliation of the observations with theoretical models, 
assuming a substantial quantity of material was accreted. 

3.3 Sodium 

The evolution of the sodium abundance as a function of lu- 
minosity for models accreting material from a 1.5 Mq com- 
panion is shown in Figure The figure also display s the 
barium-enr i ched CEMP stars from lAoki et all l|2007h and 
lAoki et al.l (|2008l ) that have measured sodium abundances. 
Most of the observed stars have sodium abundances in the 
range [Na/Fe] = — 1, though there are some extreme 
outliers which are highly enhanced in sodium, with [Na/Fe] 
around 2 or greater. 

The sodium abundance is best matched by models that 
include mixing (either by thermohaline mixing alone, or by 
the combination of thermohaline mixing and gravitational 
settling). The sodium abundance is relatively flat and shows 
no step change at or around first dredge-up. If accreted 
material remained unmixed during the main sequence, we 
would expect to see a drop of between 0.5 and 1 dex in 
[Na/Fe] at log I//Lq « 0.7. The data do not support a 
drop of this magnitude, though we caution that the data 
are a combination of a spread in metallicity and presum- 
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Figure 4. The evolution of [Na/Fe] as a function of luminosity 
for each of the model sequences. The standard models are dis- 
played with solid lines and those with only thermohaline mixing 
are displayed with dashed lines. The models with both thermoha- 
line mixing and gravitational settling are shown by dotted lines. 
From top to bottom, the panels represent the case of accreting 
0.002, 0.01 and 0.1 Mq of material respectively. In each case, ma- 
terial has been accreted from a 1.5 Mq companion and the final 
mass of the star i s 0.8 Mp^. Crosse s repre sent the Ba-rich stars of 
lAoki et ai] l|2007l ') and lAoki et al.l l|200d ). Typical random errors 
on the abundance determinations are between 0.05 and 0.15 dex 
l|Aoki et al.ll2008h . 



ably also in accreted mass and companion mass. However, 
the ejected material of the 1.5 Mq model is very sodium rich, 
with [Na/Fe] « 1.5 and any unmixed material should stand 
out. The models with thermohaline mixing bracket the bulk 
of the observations quite well, with the model where just 
0.002 M0 is accreted matching the lower observations while 
the model where 0.1 is accreted lies towards the upper 
observations. 

An alternative explanation is that the observed systems 
all have accreted material that is less sodium rich than our 



model. This would suggested that they all (excluding the 
extremely Na-rich outliers) had companions of lower mass 
than 1.5 Mq. This is possible but it would require the ma- 
jority of these binary systems to have mass ratios close to 
unity. If this is the case, it would be difficult to draw any 
conclusion about mixing from the sodium abundances. 

One may question whether the initial abundances are 
appropriate. If we look at the metal-poor but car bon-normal 
stars in the SAGA database (|Suda et al.ll2008l ). we find a 
base level of sodium to be around [Na/Fe] = 0.5, though 
there is evidence that this increases with decreasing metal- 
licity and there is some scatter. This is comparable to the 
[Na/Fe] values observed in the CEMP stars. This reinforces 
the need for mixing as the AGB material is considerably 
more sodium enhanced than this, by about one dex. 

What about those objects with enhanced sodium? The 
two most sodium-rich objects in th e sample ar e CS 2 9528- 



028 which has [Na/Fe] = 2.33 llAoki et all 12007) an d 



SDSS 1707+58 which has [Na/Fe] = 2.71 (|Aoki et al.l(2008l ) 



It is possible that these two objects had companions that 
were more massive than the 1.5 M© companions considered 
so far. In order to produce sodium during the AGB phase of 
the primary, we need first to make ^'^Ne. This is produced 
in the intershell via the reactions: 



*N(a,7)^®0(Q,7)^^Ne. 



(2) 



The addition of a proton to this ^^Ne produces ^^Na and 
this occurs in the hydrogen buring shell. The more massive 
AGB stars undergo more thermal pulses and thus can dredge 
up more ^^Ne. They also tend to have deeper convective 
envelopes, which can lead to more ^"^Na production via hot 
bottom burning. 

We have also run the model sequences accreting ma- 
terial from companions of 2 and 3.5 M©, using the ejecta 
abundances as computed by SG08 (some of which are dis- 
played in Table [l|. The ejecta from the latter model has 
[Na/Fe] = 2.50, which is greater than the observed value for 
CS 29528-028, but lower than that for SDSS 1707-^58 by 
about 0.2 dex. Thus we can reproduce the observed sodium 
abundances, but only if the accreted material is not mixed 
into the secondary. The carbon and nitrogen abudances of 
the 3.5 M0 model give [C/Fe] = 2.18 and [N/Fe] = 3.10. 
These are in reasonable agreement with the measured [C/Fe] 
and [N/Fe] values in CS 29528-028, which are 2.77 and 3.58 
respectively. The magnesium abundance is also in reasonable 
agreement as well (see below). It seems plausible that this 
object was once the companion to an AGB star of around 
3.5 Mq. We are forced to assume that the accreted material 
remained unmixed with the star, unless an order of magni- 
tude more of each element were present in the ejecta. This 
seems unlikely as discussed in Section U) 

The data points appear to show that the sodium abun- 
dance begins to rise once the luminosity reaches log L/Lq ~ 
2. Prior to this, [Na/Fe] is around 0.4 while after this point 
there is a marked increase. This trend should be treated 
with caution as there are few data points, but it may point 
to the existence of extra mixing on the red giant branch. 
This would require a non-canonical mixing mechanism that 
could reach a stellar layer deep enough for the NeNa cycle 
to be active. Proton captu res on to ^'^Ne begin at a temper- 
ature of around 2 x lO'' K (|Arnould et aLlll999l '). The ejecta 
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of the 1.5 M© model is rich in ^'^Ne and so there is enough 
raw material to produce the necessary sodium. 

Let us assume that the extra mixing that occ urs on the 
iant branch is via the ^He mechanism of lEggleton et al.l 
20061 ) and that we can reproduce th e effects using ther- 
mohal ine mixing following the work of ICharbonnel fc ZahnI 
(|2007l ). In order to do this, we need to increase the value 
of the thermohaline mi xing coefficient by about two orders 
of magnitude above the lKippenhahn et al.l l|l980l) value. We 
have re-run the post-accretion sequence of the case in which 
0.1 M0 of material is accreted from a 1.5 M© companion. 
We find that while an increase in [Na/Fe] does occur, it is of 
an insignificant amount (about 0.04 dex). In addition, there 
is a depletion of [C/Fe] at the level of about 0.1 dex and an 
increase in [N/Fe] of about 0.2 dex. The data do not support 
the occurrence of extra mixing in CEMP stars o n the upper 
part o f the giant branch (see iDenissenkov fc P insonncault 
l2008bl . for further discussion). Given the spread in the ni- 
trogen data is of the order of 1-1.5 dex, it would be easy to 
hide such a change in among the noise. 

Can we enhance this process? Enhancing the rate at 
which thermohaline mixing occurs in the giant branch is pos- 
sible, but likely to lead to the undesirable consequences that 
carbon is substantially depleted and nitrogen substantially 
enhanced. There is no observational evidence that this hap- 
pens. Also, the work of lCharbonnel fc ZahnI (|2007f ) seems to 
indicate that the efficiency of mixing we are using is cor- 
rect. So how do we produce more ^^Na? For a given tem- 
perature and density, we can increase the rate of produc- 
tion by increasing the abundance of '^'^Ne (this could be 
caused by having more efficient TDUP and/or more ther- 
mal pulses with TDUP). This element is one of the most 
abundant in the ejecta from an AGB star and it will not be- 
come severely depleted by p-burning reactions. A star which 
contains more ^'^Ne in its envelope after first dredge-up will 
produce a greater amount of sodium during extra mixing 
towards the tip of the giant branch. 

To test this we have taken one of our models after it has 
passed through first dredge-up and artificially enhanced the 
■^■^Ne abundance in the envelope by a factor of 10. Boosting 
the abundance of ^^Ne by such a large amount is somewhat 
extreme but within the variation of model predictions from 
different codes, which can vary by an order of magnitude 
(see e.g. figure 2 in SG08). The model is then evolved to the 
tip of the giant branch. We have taken the model from the 
sequence where 0.1 M0 of material from a 1.5 M0 companion 
has been accreted. Consequently the material is somewhat 
■^^Na-rich. The results are shown in Figure [S] We note that 
we can get an upturn in [Na/Fe] at the correct luminosity 
and a substantial increase in the abundance of ^^Ne in the 
accreted material does allow [Na/Fe] to a level comparable 
to that observed (the dashed line in Figure [5}. If we re- 
duce the Na abundance to a value that better represents the 
[Na/Fe] abundance after first dredge-up (i.e. we assume we 
have not accreted as much materialj, or the material is not 



* We must still satisfy the requirement [C/Fe] > +1 and this 
will indeed be met provided that the accreted mass is no less 
than about 0.01 Mq. Note that this is roughly the minimum ac- 
creted mass required for a giant to appear carbon-rich as in such 
a star the accreted material has at least been diluted during first 
dredge-up (dilution could have happened earlier, via thermoha- 
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Figure 5. The evolution of the surface [Na/Fe] abundance as a 
function of luminosity for the stellar models. The solid line dis- 
plays the model with enhanced thermohaline mixing. The dashed 
line shows the model in which the the ^^Ne abundance in the 
envelope has been increase by a factor of 10. The dotted line 
shows the model with reduced ^^Na abundance and is a better 
representation of the data. Crosses represent the Ba-rich stars of 
^Aoki et al. (2007) and Aoki et al. (2008). Typical random errors 
on the abundance determinations are between 0.05 and 0.15 dex 
jAoki et al.ll2008l) . 



as Na-rich as the AGB model predicts), we get much better 
agreement with the data. We stress that this is a simple test 
based on arbitrarily changing the Ne and Na abundances. 
Whether a stellar model (or indeed a real star) can produce 
the necessary abundances is another matter entirely. 



3.4 Magnesium 

The evolution of the surface magnesium abundance as a 
function of luminosity for models accreting 0.002, 0.01 and 
0.1 Mq of material from a 1.5 Mq companion are shown in 
Figure |6] The abundance plotted is the sum of the three sta- 
ble magnesium isotopes ^'*Mg, ^^Mg and '^^Mg. Of these, the 
AGB material is mostly enhanced in ^^Mg (with slightly less 
^''Mg) as a result of the activation of the ^^Ne(Q, n)^^Mg re- 
action in the intershell and the subsequent dredge-up of this 
material. Of the displayed models, only the models without 
thermohaline mixing pass through the observed data points 
and this only occurs while the stars are on the main se- 
quence. Once first dredge-up has occurred and the accreted 
material has become mixed, the models all skirt the bottom 
edge of the observations. 

This discrepancy is readily explained by our choice of 
solar-scaled initial abundances. In low metallicity environ- 
ments, stars appear to be enhanced in a-elements and we 
could reasonabl y expect [Mg/Fe] to be around 0.3 a t this 
metallicity fe.g. lCavrel et al.li2004l : ICohen et al.ll2004l ). This 
is presumably only due to an enhancement of the a-element 
^■^Mg. As this isotope undergoes only minor processing dur- 
ing the AGB (at least for low-mass stars like the 1.5 Mq 



line mixing). A main-sequence object whose accreted material has 
not been mixed may appear C-rich even if it has accreted far less 
material, as can be seen in figure 4 of SG08. 
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Figure 6. The evolution of [Mg/Fe] as a function of luminosity 
for each of the model sequences. The abundance is the sum of 
the three stable magnesium isotopes '^^Mg, ^^Mg and ^^Mg. The 
standard models are displayed with solid lines and those with only 
thermohaline mixing are displayed with dashed lines. The mod- 
els with both thermohaline mixing and gravitational settling are 
shown by dotted lines. From top to bottom, the panels represent 
the case of accreting 0.002, 0.01 and 0.1 Mq of material respec- 
tively. In each case, material has been accreted from a 1.5 Mq 
companion and the final mass of the star i s 0.8 WLp) . Crosse s repre - 
sent the Ba-rich stars of lAoki et al] ||2007| ') and lAoki et al.l ll2008ll . 
Typical random errors o n the abundance determinations are be- 
tween 0.05 and 0.15 dex llAoki et al.ll2008l) . 



case considered here), an increase in the initial abundance 
of ^''Mg would shift the tracks to higher [Mg/Fe] and also re- 
duce the magnitude of the dilution (irrespective of whether 
the dilution is caused by thermohaline mixing on the main 
sequence or convective mixing during first dredge-up) which 
would provide better agreement with the data. 

Figure [7] demonstrates this point. In the accreted ma- 
terial, ^^Mg is the dominant magnesium isotope, followed 
by ■^^Mg (the ratios of the magnesium isotopes in the AGB 
ejecta are ^^Mg/^'*Mg= 1.05 and ^'^Mg/^'^Mg^ 0.59.). Be- 
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Figure 7. The evolution of the magnesium isotopes as a function 
of luminosity for the model in which 0. 1 Mq of material is accreted 
from a 1.5 Mq companion and only canonical (i.e. convective) 
mixing processes are included. The tracks begin from the end 
of the accretion phase. The abundances (by mass fraction) of 
^"•Mg, -^^Mg and ^^Mg in the pristine stellar material are 2.59 X 
10-*^, 3.40 X 10-'' and 3.90 X lO"'^ respectively. Note that '^^Mg 
is the dominant isotope of Mg in the accreted material (by a 
very small margin) prior to first dredge-up, which occurs between 
logj^Q L/Lq of around 0.7 to 1.5. 



cause ■^'^Mg is little affected on the AGB, the abundance in 
the accreted material is closer to that of the pristine matter. 
When the accreted material is mixed into the secondary, the 
surface abundance of ^*Mg sufi'ers less of a depletion than 
that of the ^^Mg. The ^^Mg abundance changes by about a 
quarter, while the ^^Mg abundance is reduced by a factor 
of 4. We would expect that if we boosted the abundance of 
^■^Mg, it would dominate the total magnesium abundance 
and lead to a much reduced change in [Mg/Fe] when the 
material becomes mixed. 

In order to test this, we have re-run the case of accretion 
of 0.1 Mq of material from a 1.5 Mq companion, with ther- 
mohaline mixing and gravitational settling included. In this 
model, we have arbitrarily increased the ^*Mg abundance 
in the AGB ejecta and in the secondary by a factor of 2.5. 
This gives the secondary a magnesium-to-iron abundance 
of [Mg/Fe] — -1-0.4, which is comparable to that of carbon- 
normal EMP stars. In this model, the accreted material now 
has [Mg/Fe] = 0.68 and upon mixing, an equilibrium value 
of [Mg/Fe] = 0.4 is restored. The model now passes through 
the middle of the data, providing much better agreement, 
as shown in Figure [H] 

There is little evidence for any variation in the magne- 
sium abundance as a function of luminosity, though there 
is considerable spread in the data. This does not sit well 
with stellar models that do not include thermohaline mix- 
ing of the accreted material, as these would be expected to 
show a step-change in the Mg abundance at first dredge-up 
owing to the dilution of accreted material at this point. A 
step-change of 0.5 dex, as seen in the original solar-scaled 
standard models (i.e. those with convective mixing only), 
ought to be detectable among the spread in the data. How- 
ever, in an a-enhanced model the drop in the Mg abundance 
that occurs at first dredge-up would be reduced, as there is 
more ^''Mg in the pristine material. Consequently, the drop 
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Figure 8. The evolution of the surface [Mg/Fe] ratio as a function 
of luminosity when 0.1 Mq of material from a 1.5 Mq companion 
is accreted. Both thermohaline mixing and gravitational settling 
are included. The solid line displays the standard case while the 
dashed line is for the model with an enhanced ■^'^Mg abundance 
in both the secondary and the accreted material. 



in [Mg/Fe] at first dredge-up would become comparable to 
the scatter in the observations and therefore may not be 
evident. 

There are several objects that are extremely magnesium 
rich. The 4 most Mg-rich are: CS 29498-043, CS 29528-028, 
HE 1447+0102 and HE 0039-2635. We can reconcile the ob- 
servations of CS 29528-028 with an object that has accreted 
material from a 3.5 M© companion, but we have to assume 
the material remains unmixed on the main sequence. This 
object has [Mg/Fe] = 1.69, compared to [Mg/Fe] = 1.52 in 
the AGB ejecta. The other, more luminous, stars - which 
must have undergone first dredge-up - are difficult to ex- 
plain. This is discussed in more detail in Section |5] Here 
we simply note that a star that has passed through first 
dredge-up must necessarily have undergone some sort of ex- 
tensive dilution, even if the accreted material had remained 
unmixed on the stellar surface during the main sequence. 
This would require the accreted material of have an abun- 
dance of magnesium that is considerably higher than the 
3.5 Mq model and presumably a large amount of material 
would have to be accret ed to minimise the effe c ts of d ilution. 
This in not unfeasible: iKarakas fc Lattanziol l|2007h have a 
total magnesium abundance (which is dominated by ^^Mg) 
of around 10"'^ by mass fraction in their 3M0, Z=0.0001 
stellar model. Note that the required abundance in the ac- 
creted material increases as the mass of the accreted mate- 
rial decreases. 

It must be pointed out that even though an AGB 
model undergoes hot bottom burning, it can still be both 
a CEMP star and a nitrogen-e nhanced metal-poor star 
(NEMP). NEMPs are defined by Ijohnson et all l|2007l ') as 
having [N/Fe] > -hO.5 and [C/N] < -0.5. The star CS 29528- 
028 has [C/Fe] = 2.77 and [N/Fe] = 3.58 making it both a 
CEMP and a NEMP star. The ejecta of the 3.5 Mq model 
of SG08 has [C/Fe]= 2.18 and [N/Fe] = 3.19 which makes a 
reasonable match to CS 29528-028's abundances. 

As a caveat to the above, we note that we have assumed 
the iron abundance is [Fe/H] = -2.3, consistent with the stel- 
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Figure 9. Schematic representation of mixing after accretion. 
A star of mass Mi^it and initial composition Xq accretes Mace 
of material of composition Xacc from its companion. After first 
dredge-up, the envelope (of mass Mcnv = Mo + Macc) has become 
mixed to a homogenous compostion of Xonv 



lar models used. In fact, CS 29498-043 is considerably more 
metal-poor than this, with [Fe/H] = -3.54 (|Apki ct al. 200'^). 
This would reduce the Mg abundance required in the ejecta 
considerably. However, the other objects are of a compara- 
ble metallicity to the models presented here so this does not 
affect our conclusion on the magnesium enhancements of the 
accreted material. 



4 DISCUSSION 

We can do a simple calculation to determine the level of 
enrichment required to match the observed constraints. In 
the simplest case, we assume that there is no burning (i.e. 
the abundances are unaffected by extra mixing on the giant 
branch) of material and that material is mixed during first 
dredge-up only. The abundance in the envelope after mixing 
is determined by the depth of first dredge-up. At the deepest 
point of envelope penetration, the mass of a given isotope 
contained within the envelope is just XonvMenv where X^m 
is the mass fraction of the isotope in the envelope and Mcnv 
is the mass in the envelope. By conservation of mass, this 
must be equal to the sum of the mass of that isotope in 
the accreted matter XaccMacc plus the mass of that isotope 
contained in the original star, down to the maximum point 
of envelope penetration, XqMq (see Figure[9]for a schematic 
representation). This can be re-arranged to give: 

^onvAfcnv — XqMq 



X, 



Ma, 



which we may use to work out the required abundance of a 
given isotope in the AGB ejecta which will give the observed 
abundance on the giant branch, given an ass umed mass of 
accreted material. Stellar models (e.g. SG08, ^Lugar o et al.l 
2008) suggest that the convective envelope penetrates to a 
depth of around Mem ~ 0.45 Mq in a star of 0.8 Mq at 
Z = 10"*. 

In the case of the fiuorine-rich star HE 1305+0132 which 
has [F/Fe] = +2.9 and [Fe/H] = -2.5, the observed abun- 
dance of fiuorine implies a mass fraction of Xf — 1.65xl0~®. 
Using the above formula with a solar-scaled fiuorine abun- 
dance yields Xacc ~ 7.41 x 10~* if we assume only 0.001 Mq 
of material is accreted, while this value is one hundred 
times lower if we accrete 0.1 Mq of material. The latter case 
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is almost reconcilable with the SG08 models, which give 
Xf « 10"^ for stellar models of 1.5 and 2 Mr. The mod - 
els of iKarakas fc Lattanzid ([2007') and iLuearo et al.l (|2008l ) 
give about an order of magnitude more fluorine giving better 
agreement with the observations. 

The occurrence of thermohaline mixing does not sub- 
stantially affect the fluorine abundance required in the AGB 
ejecta. On the assumption that material mixes to the max- 
imum depth predicted by SG08 (i.e. around 0.7 M© of the 
star becomes mixed), the required mass fraction in the ac- 
creted material changes by less than a factor of two. The high 
fluorine abundance associated with only a small amount of 
accreted material is difficult to reconcile with stellar mod- 
els and it therefore seems likely that this object received 
a considerable quantity of AGB materia l. Sim ilar conclu- 
sions were also reached by iLugaro et all l|2008l) , who sug- 
gest that this object probably received between 3-11 per 
cent (0.05-0.12 Mq) of the material ejected by a 2M0 com- 
panion. These authors also point out that such extremely 
F-rich CEMP stars ought to be rare, though F-enhancement 
likely occurs alongside C- and s-processes enhancements and 
should be sought out as a means of confirming the AGB mass 
transfer scenario. 

Repeating this analysis for the highest observed sodium 
and magnesium enrichments, we find required abundances 
vary from XNa = 3x10"'' -10"^ and Xms = 8x lO'^-lO"^. 
Again, the high abundances are irreconcilable with stellar 
models. It seems highly unlikely that stellar models would 
produce sodium or magnesium at mass fractions of around 
1 per cent! Current stellar models predict the upper limit 
on the abundance s of these two elements to be in the range 
10"" - 10"^ (e.g. lHerwiell2004l : lKarakas fc Lattanzid l2007l . 
and SG08), depending on the input physics and the mass of 
the stellar model. This upper limit is roughly consistent with 
our simple calculations for the required sodium and magne- 
sium abundances in the accreted material. It is tempting 
to suggest that many CEMP stars do not accrete a great 
deal of material from their companions, on the grounds that 
we do not see many highly enriched objects (regardless of 
whether that enrichment is in terms of F, Na or Mg). 

The models computed in this work have not included 
the effects of radiative levitation. This would be expected 
to temporarily raise the surface abundance of some of the 
heavy elements, p articularly Na and Mg (see figure 4 in 
iRichard et"al]|2002i '). It would also raise the Fe abundance, so 
the exact effect on the ratio [X/Fe] is not straightforward to 
predict. Radiative levitation will not change the qualitative 
picture presented herein because it is only effective over long 
(i.e. Gyr) timescales just like gravitational settling. It does 
not substantially afi'ect the settling of helium and hence the 
/^-barrier that may develop. As it is this that determines the 
extent to which thermohaline mixing can occur and because 
thermohaline mixing operates over a short timescale relative 
to the settling processes, the qualitative picture presented 
herein should not change. In addition, if an extra turbulent 
mixing process is needed to explain the differences in the 
iron and magnesium abundances that exist betwee n tur n-off 
and evolved stars in globular clusters (iLind et al.ii2008l ), the 
omission of radiative levitation will have less of an effect as 
extra mixing reduces the impact of settling and levitation. 
The predictions for the giant branch will be unchanged be- 
cause the material in the envelope is thoroughly mixed by 



convection; it is only the main-sequence predictions that will 
be affected. 

A wealth of information is also available for several 
heavy species not included in our network. Several authors 
have compared the s-process element abundance trends to 
predictions to stellar models and come up with factors re- 
flecting the d egree of dilution that must have occurred i n 
the star (e.g. iThompson et aLllioOsI : iRoederer et al.ll2008l ). 
In particular, comparing these heavy metal abundance pat- 
terns in turn-off stars to model predictions may help to elu- 
cidate the nature and extent of any mixing on the main- 
sequence. We intend to pursue this avenue in future work. 

The work discussed herein only applies to those CEMP 
stars that have formed via the accretion of material from an 
AGB companion. This likely applies to the CEMP-s stars, 
which account for the majority of the CEMP p opulation 
(around 70 per cent, according to lAoki et al.|[200a ). The re- 
maining CEMP stars are likely to have formed via a different 
mechanism. Other proposed formation scenarios invoke the 
ejecta of rapidly rotating ma ssive stars (|Mevnet et al.ll2006l ) 
or faint type II supernovae (jUmeda fc Nomotd 20051 ). It is 
also possible for extremely metal-poor stars to enrich them- 
selves through a dual core flash, where protons are ingested 
into the convective region that is driven during the core 
He-fl ash at the tip of the giant branch (e.g. iFuiimoto et al.l 
'l99Cf). This o nly happens for stars with metallicities below 
[Fe/H] ^ — 5 IjCampbell fc Lattanzio 2008 ). which is consid- 
erably more metal-poor than the objects discussed in this 
work. 



5 CONCLUSION 

The measurement of a high Li abundance in an unevolved 
CEMP star precludes the possibility of extensive mixing of 
accreted material. This could either be because there is no 
thermohaline mixing or because only a small quantity of 
material is accreted (likely a few thousandths of a solar 
mass). The inclusion of gravitational settling does not sub- 
stantially increase the very stringent limit imposed on the 
amount of mass that can be accreted. It is possible that 
the Li abundance in accreted AGB material is substantially 
higher than computed at present owing to the occurrence 
of some additional mixing mechanism during the thermally 
pulsing asymptotic giant branch. However, the Li enhance- 
ments required to make the thermohaline mixing models fit 
the observations would have to be at least two orders of 
magnitude above the current yield. This seems an unlikely 
level of enhancement. 

Comparing measured sodium abundances to the models 
suggests that the majority, but not all, stars undergo mixing 
during the main sequence. The abundance in the AGB ejecta 
is significantly higher than that observed in most CEMP 
stars. Unless sodium is produced to a much lower degree 
in AGB stars than the models predict (or the stars we are 
observing all had companions of much lower mass than we 
have assumed), it is difficult to reconcile models without 
main-sequence mixing with the observed abundances. How- 
ever, it is difficult to explain the extremely sodium-rich ob- 
jects with models that do include thermohaline mixing as 
this would require the accreted material to have a sodium 
abundance considerably higher than current model predic- 
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tions. It therefore seems prudent to argue that some stars 
mix their accreted material while others do not or that there 
is a spread in the efficiency of the mixing of accreted ma- 
terial. The cause for this is unclear and merits investigat- 
ing. One possible suggestion is that rot ation may play a 
role. IPenissenkov fc Pinsonneauld ()2008ah have pointed out 
that rotationally-driven horizontal turbulence may suppress 
thermohaline mixing. If the secondary stars are rotating at 
different rates (or the angular momentum content of the 
accreted material varies from star-to-star) then some stars 
may experience thermohaline mixing whilst others do not. 

The flatness of the observed magnesium abundances 
also does not support models without thermohaline mixing 
on the main sequence. However, this may be because the 
models used are not a-enhanced and consequently are likely 
to be less rich in ^*Mg than the observed stars. AGB and 
secondary models with a-enhanced compositions and abun- 
dances should certainly be investigated. As the magnesium 
from AGB stars is predominantly ^^'^®Mg, additional ^^Mg 
in the stellar material could easily mask the dilution of the 
AGB material and hence we would not be able to see any 
changes due to mixing events. However, if the CEMP stars 
have received a substantial quantity of magnesium from an 
AGB companion, it is likely to be in the form of the neutron- 
rich isotopes ^^Mg and ^^Mg. Hence the signature of accre- 
tion from such a companion may be detectable in the ratios 
of ^^'^®Mg to ^^Mg. Such measurements have been made 
in metal-poo r but carbon-normal stars but the results are 
inconclusive (lYong et al.ll2003l : iMelendez fc Cohenll2007l '). 

The cases of extreme enhancement of F, Na and Mg on 
the giant branch argue against the accretion of only small 
quantities of material (i.e. thousandths of a solar mass) in 
these cases. The required abundances in the AGB ejecta 
would be beyond the current range of model predictions. 
The accretion of much greater amounts of material would 
be expected to show more extensive mixing on the main 
sequence if thermohaline mixing is efficient. 

No clear picture arises from the elements studied here. 
The detection of lithium in turn-off objects seems to rule out 
extensive mixing for these cases. Some highly sodium and 
magnesium enriched turn-off objects can only be explained 
if accreted material remains unmixed during the main se- 
quence. However, measurements of [Na/Fe] seems to sug- 
gest that most CEMP stars do efficiently mix their accreted 
material. A more comprehensive study involving a greater 
number of elements may help to elucidate what processes go 
on in these objects. 
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